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ABSTRACT 

We present Spitzer IRS observations of a complete sample of 57 SDSS type-1 quasars at z ~ 1. 
Aromatic features at 6.2 and/or 7.7 are detected in about half of the sample and show profiles 
similar to those seen in normal galaxies at both low- and high-redshift, indicating a star- formation 
origin for the features. Based on the ratio of aromatic to star-formation IR (SFIR) luminosities for 
normal star- forming galaxies at z *~ 1, we have constructed the SFIR luminosity function (LF) of z 
~ 1 quasars. As we found earlier for low-redshift PG quasars, these z ~ 1 quasars show a flatter 
SFIR LF than do z ~ 1 field galaxies, implying the quasar host galaxy population has on average 
a higher SFR than the field galaxies do. As measured from their SFIR LF, individual quasar hosts 
have on average LIRG-level SFRs, which mainly arise in the circumnuclear regions. By comparing 
with similar measurements of low-redshift PG quasars, we find that the comoving SFIR luminosity 
density in quasar hosts shows a much larger increase with redshift than that in field galaxies. The 
behavior is consistent with pure density evolution since the average SFR and the average SFR/BH- 
accretion-rate in quasar hosts show little evolution with redshift. For individual quasars, we have 
found a correlation between the aromatic-based SFR and the luminosity of the nuclear radiation, 
consistent with predictions of some theoretical models. We propose that type 1 quasars reside in a 
distinct galaxy population that shows elliptical morphology but that harbors a significant fraction of 
intermediate-age stars and is experiencing intense circumnuclear star formation. 
Subject headings: galaxies: active — galaxies: nuclei — galaxies: starburst 



1. INTRODUCTION 

Since its discovery, the relationship between the 
black-hole (BH) mass and galaxy bulge properties 
(IKormendv fc Richston el 119951: [Magorrian et all 119981 : 
iGebhardt et all 120001: " iFerrarese fc MerrittJ l2000f ) has 
played a critical role in our ideas about galaxy evolu- 
tion and BH growth. Active galactic nuclei (AGNs), the 
manifestation of accretion onto BHs, are the main sites 
of BH growth. The demography of local galaxies sug- 
gests that most - perha ps all - massive galaxies hos t 
BHs at their centers (e.g. IKormendv fc Richstond ll995L 
The good match between the local BH mass density and 
the accreted BH mass density in AGNs suggests that the 
AGNis an indispensable phas e of ga l axy evolution (e.g. 
Soltanl [l98l lYu fc Tremaind 120021: lAller fc Richstonel 
2002t IShankar et al.ll2004j : iMarconi et al.ll2004D . The in- 



terplay between star formation and AGN activity may 
play a key role in the establishment of the BH-bulge cor- 
relation. However, in quasars (luminous AGNs), obser- 
vational constrains on the BH/star- formation relation- 
ship are rare due to the difficulty in measuring the star 
formation rate (SFR) around the bright quasar nuclei. 

The host galaxy morphology and colors provide the 
first insight into the level of star formation. At low 
redshift (z < 0.5), decades of effort have been devoted 
to imaging the host galaxies o f UV 7 optically-se lected 
type 1 quasars (IHutchings et allll984tlSmith et al.lll986l : 



McLeod fc Riekdll995HBahcall et alJll997t Dunlop et all 
2003t iFlovd et all 120041: IGuvon et all 12009) . Important 



technical developments, especially with space-based and 
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adaptive optics (AO) instruments, have provided obser- 
vations with dramatically improved spatial resolution, 
stable PSFs and large dynamic range, which are keys 
to extracting quantitative information on host galax- 
ies surrounding the bright AGN. Radio-loud quasars are 
found to be associated with luminous early-type galaxies 
with structure paramet ers indistinguishable from those 
of qu iescent ellipticals (|Smith et all 119861 : iDunlop et all 
I2003T ) . Radio-quiet quasars reside on average in less lu- 
minous galaxies with diverse morphologies that show 
some dependence on the host and nuclear luminosity, 
i.e., more frequent ly of early type for the more luminous 
nucle i and hosts jlMcLeod fc Riekd 119951: iBahcall et all 
[19971 IFlovd et all 120041 : IGuvon et all 12009)" At M B 
< -23, roughly 50% of quasar hosts show normal el- 
liptic al morphologies (|Bahcall et all Il997t IGuvon et all 
12009) . The fraction of bright qu asar hosts with on- 
going merging feat ures is 10-30% (Ba hcall et all 119971 : 
IGuvon et all [2006). However, this merger fraction 
may be under-estimated especially at the late merg- 
ing stage where the tidal tails become relatively faint. 
Such a bias has been demonstrated by the discov- 
eries of merger remnants in quiescent e lliptical hosts 
([Canalizo et all 120071 : iBennert et al.ll2008D. The merger 
fraction for near-IR-selccted quasars dCutri et all 2002) 
and far-IR excess quasars (Canalizo & Stoc kton] 12001 ) 
is much higher (>50%) (Canaliz o fc Stockton! 12001 : 
IHutchings et al.ll2003l : iMarble et al.ll2003l ). However, it is 
unclear if this is because tidal features are relatively eas- 
ily detected in IR-selected quasars with low nuclear/host- 
light contrast or if there is an evolutionary difference 
between quasars selected in the UV/optical and those 
selected in the IR. 

The study of stellar populations in quasar hosts using 
UV/optical spectra has provided another significant 
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constraint on their star formation properties, mainly 
through detecting the presence of the young- (<100 
Myr), intermediate- (~0.1-lGyr) or old- (~10Gyr) 
stellar populations. However, in both type-1 and type-2 
quasars, the nuclear radiation strongly dilutes the 
young stellar features, such as the PCygni UV lines 
and the 4650A WR bump, and dominates the line 
emission of SFR tracers used for normal galaxies. The 
mean-age-indicator 4000A break and intermediate-age- 
indicator (~0.1-1 Gyr) H5 absorption line also suffer 
from dilution by nuclear light. Despite these difficulties, 
and inconsistent with the implication of the early-type 
host morphologies, an intermediate-age stellar pop- 
ulation is present at a significant level. Off-nuclear 
(>~15 kpc) UV/optical spectra of quasars indicate that 
old stellar popul ations dominate ijHughes et al.l 120001 : 
iNolan et al.ll200ll ). However, with much deeper spectra 
a significant fraction (10% in mass f raction) of an 
intermediate-age population is rev ealed (|Canalizo et al.l 
[2fM lOanalizo fc Stocktonl Hp|). In type-2 AGNs 
and less luminous type-1 objects where the nuclear 
radiation contamination is less severe, the presence of an 
intermediat e-age stellar populat i on ha s been extensively 
confi rme d dKotilainen fc Wardl 1 1 9941 : iRonnback et all 
19961: iBrotherton et al.l )1999t iKauffmaim et"afl 
200a iJahnke et all 12004 (Vanden Berk et all 1200a 
Jahnk e~et al.l 12007). Furthermore, it seems there is a 



decrease in the mean stellar age w i th increasing AGN 
lumin osity ((Kauffmann et all 120031 : IVanden Berk et all 
I2006T ). 

Techniques that are more suitable for detecting ongo- 
ing/recent (< 100 Myr) star formation have been em- 
ployed recently, although their accuracy is still not com- 
parable to those for normal galaxies. As the mate- 
rial reservoir of star formation, the cold molecular gas 
mass roughly correlates with the level of star forma- 
tion. It has been fo und that quasar h o sts are rich 
in cold molecular gas (IScoville et all 120031 : 1 Evans et al.l 
120061 : iBertram et al.ll2007D . implying significant ongoing 
star formation activity. With space-based and AO in- 
struments, spatially-resolved images of SFR tracers are 
promising ways to quantify the SFR and characterize the 
spatial distribution of the star-forming regions. For ex- 
ample, the extended Paa emission in PG1126-041 is most 
likely associated with an intense nuclear (lOOpc) s tar- 
burst embedded in the old bulge (|Cresci et al.ll2004l ). 

A new method to measure the SFR of local quasar 
hosts is to use the mid-infrared aromatic bands 
(IQgle et al.H200a iSchweitzer et al J 120061 1 Shi et al.ll2007l: 
iFu fc Stocktonl l2009f ) measured with the Spitzer Space 
Telescope. In field galaxies, these emission bands are 
universally bright in roughly solar metallicity regions 
of star formati on, with well understood relationships to 



the SFR (e.g. iRoussel et~all l2001|; pale fc Helou] 



Wu et alJ 120051: iBrandl et al.l 120061 : ISmith et al.1 



2002; 



20071 : 



Engelbracht et al.l 120081 ). Their carriers appear to be 
destroyed by the harsh radiation f i elds around high - 
luminosity AGNs (jGenzel et al.lll998tlS~p"oon et al.ll2007h 
and thus any detected aromatic emission is likely from 
the host galaxy 

In IShi et all (|2007h , we provided further evidence for 
the star-formation origin of the aromatic features in 
quasar hosts: (1) both the overall shape of the aromatic 
features and the distribution of the 7.7/im/11.3^m fea- 



ture ratio are similar to those in star forming galaxies; 
(2) the larger the equivalent widths of the aromatic fea- 
tures, the stronger the far-IR emission in the global SEDs 
of quasars; (3) quasars lie on the trend of normal galaxies 
in the plane of the molecular gas mass and the aromatic- 
derived SF IR luminosity. By measuring the aromatic 
feature fluxes in a large quasar sample (^200 objects in- 
cluding PG, 2MASS & 3CR), we obtained a quantitative 
census of star formati on activity in low-redshift quasar 
hosts (|Shi et al.l l2007f) . We showed that quasar hosts 
have a flatter star formation infrared (SFIR) luminos- 
ity function (LF) than field galaxies, i.e., the quasars lie 
in very luminous star forming galaxies more often than 
would be the case for a random sampling of galaxies with- 
out nuclear activity. 

As a summary, type-1 quasars at low redshift (z < 
0.5) mainly reside in luminous early- type host galaxies 
but harbor a significant (>10%) fraction of intermediate- 
age stellar populations and are experiencing intense star 
formation activity. These properties are rare for nor- 
mal galaxies, implying that the quasar appears at a spe- 
cial stage of galaxy evolution. In this paper, we probe 
whether type-1 quasar host galaxies at z ~ 1 have similar 
characteristics. 

As demonstrated in our low-redshift work (|Shi et al.l 
l2007f l. the aromatic feature observed with Spitzer pro- 
vides efficient measurements of the SFR in quasar hosts 
with significantly improved accuracy. To constrain the 
star- formation/quasar interplay at high redshift and to 
characterize the cosmic evolution of star formation in 
quasar hosts, we have carried out IRS observations of 
a complete optically-selected type-1 quasar sample (57 
objects) at z ~ 1 to constrain their SFRs. The sample 
selection and data reduction are presented in § [2] The 
measurement of the aromatic feature flux is discussed in 
§ [3] We show the results in § [4] and discuss them in § [5j 
The conclusion is presented in § [6] Throughout the pa- 
per, we adopt a cosmology with H o =70 km s _1 Mpc -1 , 
fi m =0.3 and £!a=0.7. All magnitudes are defined in the 
AB magnitude system. 

2. SAMPLE AND DATA REDUCTION 

Our quasar sample (Table [TJ is a subset of the 
nearly homogeneous parent sample from Ric hards et al.l 
(2006). The completeness for this parent sample has 
been well quantified and it has been used to construct 
the quasar LF. To observe objects in low IR back- 
gro und, we selected all ta rgets at 0.8 < z < 1.0 from 
the iRichards et al.l (|2006f ) sample in three fields cen- 
tered at but significantly larger than three Spitzer deep 
fields, including the Spitzer Wide-Area Infrared Extra- 
galactic Survey (SWIRE) ELAIS-N1, SWIRE ELAIS- 
N2 and SWIRE Lockman Hole. Our three fields have 
sizes of 15.85° x9.01° centered at 16hllm00s +55°00"00', 
19.22° x 12.0° centered at 16h36m48s 41°01"45', and 
3.77° x 3.0° centered at 10h45m00s +58°00'00", respec- 
tively. Our final sample is composed of 57 SDSS quasars 
at 0.8 < z < 1.0 and is complete down to SDSS i=19.0. 
As the SDSS quasar coverage does not fill the above 
three fields, the effective area is calculated by subdi- 
viding the field into small rectangles (about 1 square 
degree) and then summing those falling within the ac- 
tual survey region. This gives an effective area of 88.0 
square degree, which is almost the same as the value (83.0 
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TABLE 1 

Aromatic Features in z ~1 Quasar Hosts 



Sources 


z 


mi 


M g 


^6.2PAH 
[10~ 15 crg s -1 cm -2 ] 


^Y.7PAH 
[10- 15 erg s" 1 cm" 2 ] 


r SF 
TIR 

[1O U L ] 


^5-6//m 

[10 n L o ] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


blJbblOooo5.oo+5/5ol4. 7 


0.96 


18.56 


-24.87 


< 2.0 




< 3.07 


1.19 


bJJbbl04Zoy.00+5ooZol.0 


1 nn 
1.00 


iy.00 


O/l A A 


<C z.Z 




< 4.Z4 


n a a 
U.44 


cmQCinf;nnn 01 i r;cinn/i o 
blJbb 100000. Zl+DoiyU4.Z 


n qq 
O.00 


1 1 Tfi 
1 ( . (0 


or; qq 
-ZO.00 


< l.O 




1 1 n 

< i.iy 


1 nn 
1.00 


qticci K/i K/io 7o i r;nr;'7'r;n o 
blJbbl5454z. /5+505/5y.z 


0.94 


19.04 


-24.35 


< 1.4 




< 1.67 


0.52 


o T"\ o o i rr tr o rr q /i i a oc o 1 t t 

bJJbbl5555o.o4+4Z55l7.7 


0.87 


18.50 


-24.75 


1.5+0.5 




1.43+ 1.17 


0.51 


bDbblooyo4.o 7 +0o0 (41 .4 


0.96 


18.39 


-25.03 


< 1.6 




< 2.25 


0.54 


bDbb Io001o\44+ooo40o.z 


0.95 


19.02 


-24.39 


3.6+0.5 




t 1 I /i 7n 

7.1z± 4.79 


0.78 


QHQQ1 o cn i /i/ic:Q/in n 
bJJbbloOolo.oO+44oo4y.O 


0.98 


19.05 


O/l A O 

-Z4.4Z 


< 2.0 




< 3.45 


0.71 


nriQC 1 ftn /I /1 1 oc i /I ft Q'zcc; o 

bJJbblo0441.oo+4oo (5b. o 


0.94 


19.09 


O/i on 
-Z4.00 


< 1.3 




< 1.56 


0.40 


CFlQG 1 ftnr;Q'7 K/1 1 Q700D1 O 

bJJbb loooz / .o4+o i zsoi.y 


n qo 
O.oZ 


1 T 

1 ( .O ( 


or; co 
-ZO.Oo 


^ 1 Q 
< l.O 




^ -\ AI 
< 1.4/ 


1 10 

1.1Z 


nr\QG i ftnft r;n r;n i at Am c\ q 

bJJbbio0ooy.oy+4i4yiy.o 


0.98 


18.97 


-24.49 


< 1.4 




< 2.10 


n C7 

0.57 


bDbbl60733. 26+451610. 6 


0.94 


18.93 


-24.47 


2.9+0.5 




5.04+ 0.49 


0.54 


bDbblo07o5.49+4b575o.o 


0.97 


19.06 


-24.38 


< 1.4 




< 1.94 


0.97 


bDbblo0ozZ.0z+415o5z. / 


0.84 


17.42 


or; "7 a 
-Z0. 74: 


< 1.4 


1 n n_i_o n 

io.y±z.o 


no_L 01c 
z.yz± Z.10 


1.78 


Qr\cc i ftnQ err: a Q i /IQrrocn O 
blJbbl00o05.40+4o5Z5y.Z 


n nQ 

o.yo 


1 / .ZO 


OC\ 1 Q 

-Z0. lo 


< 1. 7 




/ rt rip- 
< Z.ZO 


A OQ 
4.ZO 


nr\QG 1 aC\C\(\A d. A 1 /I nc C d *7 r; 

bJJbblo0y04.o4+40oo4/ .5 


0.89 


18.70 


-24.58 


3.9+0.4 


1 a n_L a n 
lo.0±4.0 


a 1 Q _L A 01 
0.1o± 4.Z1 


0.89 


bDbb 161055. 65+380305. 7 


0.83 


18.54 


-24.60 


4.1+0.4 


15.5+1.8 


c or* 1 c q 
0.Z0+ 0.00 


0.66 


bJJbbloll04.75+4bz05z.o 


0.83 


18.80 


-24.33 


< 1.4 




< 1.06 


1.82 


nriQC i a 1 i ^ A n 1 i 97nQ9n v 
bDbb lol 144. Ol+o (OooO. 7 


0.81 


17.77 


-25.31 


2.7+0.6 


11.2+1.7 


O /;r 1 -1 nn 

z.oo± i.yy 


1.15 


bDbb lol Zoz.50+40oOoz.O 


0.82 


17.16 


-25.96 


6.2+0.5 


OQ Q_l_1 Q 
Z5.0±1.0 


y.lo± 0.0/ 


3.00 


biJbbioio4z.y / +oyo roz.e 


0.98 


18.25 


-25.21 


< 1.5 




< 2.30 


1.07 


nriQC 1 ft 1 Q r; 1 Q/l l Q7y|OC;0 *7 

bDbbloloOl.o4+o/4Z5o. / 


0.81 


18.87 


-24.21 


< 1.6 




< 1.23 


0.42 


cr>QC i fti /I t;c r;ft i AP.C\1AA C\ 
bDbb 10 1400. 00+4b0 / 44. U 


n q c; 
O.oO 


is. ( y 


O/l QQ 

-Z4.oo 


1 n re 
Z.Z±0.0 


OI OJ-O Q 

Zl.ZztZ.O 


Z.Zozt 1./0 


1.14 


encci ^1 /JQ7 i c i onnQCft q 
bDbb lo loo /.Io+oy0o5b.o 


0.81 


17.89 


-25.19 


3.2+0.8 


14.0+1.7 


Q A O _L O A O 

o.4Z± Z.4o 


1.72 


CnCG1C17Ki! 7"7 I /IOOnO/1 1 

bDbblol7oo.77+4Zoyz4.1 


0.97 


18.62 


-24.84 


2.3+0.5 




A no 1 on 


0.80 


Qnocicionc qi i /ioocqo 1 
bDbbloloOo.ol+4zzooz.l 


0.93 


1 ( .94 


OC A A 

-Z0.44 


< 1.6 




< 2.01 


2.40 


bL)bbloly4o.zy+4ooyio.4 


0.85 


18.61 


-24.58 


2.0+0.6 




n 1 _l_ 1 r;T 
Z.01± 1.57 


0.40 


bDbbloz0o5.bo+4z0 ( 4z.o 


0.96 


19.08 


O/l Q£i 

-Z4.00 


< 1.9 




< 2.93 


0.48 


Gncci^ont;i 17 1 aoiaah q 
bDbb lozOo 1.1 7 +4zo44y.o 


0.99 


18.88 


O/i a i 
-Z4.01 


< 1.7 




< z.7y 


1.69 


oncci /;oniro no i zion/iO/i cr 

bDbbloz058.08+4204z4.5 


0.96 


18.84 


-24.60 


< 1.9 




< 2.93 


0.39 


bDbblozll0.oo+oblo5o.o 


0.83 


17.88 


-25.26 


< 1.3 




< 0.99 


0.79 


nr\cc OO Q£? i /I Q CO on A 

bDbblozlZo.oo+4Z5zzy.4 


0.98 


18.42 


-25.04 


3.8+0.5 




q a A _i_ r; a. i 
o.44± o.ol 


1.27 


bDbblozlo5.b5+oy545z.l 


0.98 


19.05 


-24.41 


< 1.8 




< 2.87 


0.50 


nr\C C 1 eoOO /I ftft 1 971 QHO n 

bDbblozzz4.bo+o/ loOU.y 


0.82 


19.05 


-24.05 


2.8+0.4 


10.6+1.7 


O Q Q _1_ O 1 Q 

Z.50± Z.lo 


0.48 


encci ftoo/iQ qo i /inQnoo a 
bDbb!0ZZ4o.oZ+4Uo0Zy .4 


n qi 
O.ol 


1 q nQ 

lo.yo 


O/l 1 A 
-Z4. 14 


< 1.0 




< l.zo 


n ft/i 
0.04 


o T^i c o i o o i o o n i /innocro T 

bDbb 162318. 89+402258. 7 


0.91 


17.78 


-25.55 


< 1.4 




< 1.50 


1.20 


bDbb 162330. 53+355933. 2 


0.87 


18.46 


-24.76 


< 1.5 


11.8+3.1 


Q TO 1 O / ' "7 

3.72+ 2.67 


1.15 


bDbb 162349. 46+364721. 7 


0.92 


18.49 


-24.87 


< 1.6 




< 1.93 


0.72 


bDbb!62449. 96+35085 r .9 


0.88 


18.94 


-24.32 


< 1.6 




< 1.69 


0.68 


bDbb 162658. 11+353030.1 


0.85 


18.76 


-24.43 


< 1.5 


7.7+2.1 


1 /i _i_ 1 c 
1.84± 1.46 


0.55 


bDbb!62902.59+3 ( 2430.8 


0.93 


19.03 


-24.34 


3.2+0.5 




5.60± 3.84 


0.76 


bDbb!63018. ( l+o ( iyU4.5 


n n'T 
O.y 1 


18.39 


or; nc 

-zo.Oo 


< 1.5 




< 2.04 


1.84 


SDSS163225. 56+411852. 4 


0.91 


18.42 


-24.91 


1.7+0.4 




2.07+ 1.61 


0.83 


SDSS163537.99+365936.6 


0.90 


18.35 


-24.97 


< 1.4 




< 1.53 


1.26 


SDSS163624.98+361457.9 


0.91 


18.93 


-24.40 


< 1.5 


17.8+4.5 


7.92+ 5.29 


0.33 


SDSS163656.84+364340.4 


0.85 


18.54 


-24.65 


< 1.5 


7.4+2.4 


1.72+ 1.38 


1.19 


SDSS163726.88+404432.9 


0.86 


18.52 


-24.68 


4.1+0.4 


7.9+2.3 


6.03+ 4.11 


0.66 


SDSS163926. 19+390821.4 


0.88 


18.88 


-24.37 


< 1.4 




< 1.30 


0.21 


SDSS164334.80+390102.5 


0.93 


18.91 


-24.45 


< 1.6 




< 1.97 


0.54 


SDSS164346.66+383812.9 


0.89 


18.94 


-24.35 


5.6+0.4 


13.6+3.6 


10.50+ 6.92 


0.36 


SDSS164455.12+384936.3 


0.86 


18.78 


-24.43 


< 1.5 


13.5+3.1 


4.34+ 3.06 


1.23 


SDSS164508.74+374057.0 


0.93 


19.02 


-24.35 


1.9+0.5 




2.58+ 1.94 


0.73 


SDSS164600.09+381833.3 


0.99 


18.18 


-25.30 


< 1.5 




< 2.36 


2.20 


SDSS164730.89+360101.6 


0.85 


18.78 


-24.40 


< 1.4 


7.1+2.3 


1.64+ 1.32 


0.63 


SDSS164745.05+355732.9 


0.94 


19.07 


-24.33 


1.4+0.5 




1.81+ 1.44 


0.43 


SDSS164925.07+373015.4 


0.95 


18.97 


-24.44 


4.9+0.5 




10.96+ 7.21 


0.53 


SDSS165614.00+351014.5 


0.81 


17.44 


-25.63 


3.8+0.6 


29.9+1.7 


4.37+ 3.08 


1.73 


LogL 5 _ 6Mm e [ 8.0, 10.8] 












0.18tg;" 


43+O.14 


LogL 5 _ 6Mm 6 [ 10.8, 12.5] 












n or; + l-22 


1 04+2.99 



Note. — Col.(l): Source name. Col. (2): Rcdshift. Col. (3): the i-band magnitude. Col. (4): the rest-frame g-band magnitude. 
Col. (5): the 6.2 fim aromatic feature flux and 3cr uppcrlimits. Col. (6): the 7.7 fim aromatic feature flux. Col. (7): the PAH-dcrivcd 
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Fig. 1. — Continued 



square degree) derived from the total effective area of the 
SDSS-DR3 survey and the relative number of SDSS-DR3 
quasars at 0.8 < z < 1.0 included in our sample. Al- 
though we observed quasars in low IR background, there 
are no selection criteria based on the IR fluxes of indi- 
vidual objects. 

We obtained IRS spectra of the members of this sample 
(Pi-George Rieke, PID-50196). For objects in the red- 
shift range between 0.8 and 1.0, IRS SL1 (7.4-14.5/xm) 
covers the 6.2 /jm aromatic feature and the spectral range 
required for continuum subtraction. For objects at z < 
0.9, the blue wing of the 7.7 /im feature will also be ob- 
served if it is present. The IRS staring mode was used for 
the observations and the exposure time for each object 
was 240x2 sees. 

The data reduction basically followed the IRS data 



handbook. Briefly, each BCD image was first cleaned 
using IRSCLEAN. The images observed on the same slit 
position were sigma-clipped, averaged, and used to sub- 
tract the sky for each image observed on another slit 
position. The spectra were extracted from each sky- 
subtracted BCD image with the optimal extraction al- 
gorithm using SPICE 3 and then sigma-clipped and av- 
eraged to produce the final combined spectra. We indi- 
vidually inspected each flux profile along the slit to make 
sure the extraction aperture was centered on the bright- 
ness peak. 

3. MEASUREMENTS OF AROMATIC FEATURES 

Because we can work at shorter rest wavelengths, the 
extraction of the aromatic features is relatively straight- 

3 http://ssc.spitzer.caltech.edu/postbcd/spicc.html 
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forwar d compared to that for local AGNs in iShi et all 
(|2007f ). In that case, careful subtraction of silicate fea- 
tures was needed to measure the 7.7 and 11.3 /im aro- 
matic bands and Monte-Carlo simulations were carried 
out to characterize the measurement uncertainties. In 
the current situation, for each spectrum the continuum 
is fitted with a power law shape to the flux density in the 
spectral range of 5.8-6.0 /im and 6.8-7.1 /im. We then 
fitted the continuum-subtracted spectra with two Drudc 
profiles, one centered at 6.22 /im with a FWHM of 0.186 
/um and another centered at 7.7 /im with a FWHM of 
0.53 Jim for the 6.2 and 7.7 /im features, respectively 
fsee lSmith et al.ll2007l ). For several objects, the central 
wavelength of the 6.2 /im feature is adjusted between 6.2 
and 6.3 /im to achieve the best fit. 
The feature flux was obtained by integrating the fit- 



ted profile from 1 to 20 /im. For the 7.7 /im feature, 
the flux was only measured if the rest-frame spectrum 
wavelength extends longer than 7.4 /im. The noise of 
the 6.2 /tm feature was defined within the spectral range 
between 6.0 /im and 6.4 /tm. For the 7.7 /tm feature, the 
S/N is defined in the spectral range from 7.0 /tm to the 
longest wavelength (usually < 8/tm). Aromatic features 
with S/Ns above 3 were considered to be solid detections. 
We inspected each detected feature and excluded several 
detections due to large noise around the feature. The 
final detection rate is ~37% for the 6.2 /im feature and 
~50% for any feature at either 6.2 or 7.7 /tm. The spec- 
tra and fitted features are shown in Figure [1] The fluxes 
or 3tr upper- limits are listed in Table [TJ 

For the 7.7 /im feature, a te st based on t h e four nearby 
galaxy composite spectra of ISmith et al.l (|2007[ ) shows 
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Fig. 1. — Continued 

that a fit only using the blue wing does not introduce a 
large uncertainty. On the other hand, the silicate emis- 
sion feature and sometimes [NeV |7.65 /im line c an be 
present in these type 1 quasars (jShi et all 12006). To 
quantify the effect of these two features on the 7.7 fim 
aromatic fluxes, we have used the above method to mea- 
sure the feature fluxes in all PG quasar s. They are com - 
pared to more accurate ones obtained in lShi et al.l (|2007l ) 
by deriving (/ 7 .7 M mPAH - /7%mPAH)//7.7 M niPAH> where 
/7.7/jmPAH is the flux based on the method of this pa- 
per and fj 7/imP AH i s * ne one derived in lShi et ail (|2007f ). 
This ratio has an average value of -0.2 and a 68% confi- 
dence range of -0.55 to 0.15. The differences are mainly 
caused by the deviation of the IR spectrum at A > 8 /xm 
from the interpolation based on the part below 8 /im. 
Such a deviation is due to either a change in the IR con- 



tinuum or the presence of the silicate emission feature. 
Ten objects out of fifty-seven have both 6.2 and 7.7 /im 
features and seven more have detections only of the 7.7 
/im feature. The median value of the /7.7PAH//6.2PAH ra- 
tio is 4.1±2.4. This ratio is used to obtain the equivalent 
6.2 /im feature flux for those objects with only 7.7 /im 
fluxes. 

4. RESULTS 

4.1. IR Spectra of z~l Quasars And Comparisons To 
Other Samples 

IShi et all (|2007l ) found that the profiles and relative 
strengths of the aromatic features in local AGNs are 
quite similar to those of local star-forming galaxies, im- 
plying star formation excitation of the aromatic fea- 
tures in the AGNs. Figure [2] compares the continuum- 



Shi et al. 



E 



0.3 
0.2 
0.1 
0.0 
0.3 
0.2 
0.1 
0.0 
0.3 
0.2 
0.1 
0.0 



(a) 



z = Gol. (Smith07) 
SDSS 




z=1.7 ULIRGs(Forroh08) 
SDSS 




(e) 



'1 1. 1 1 1 ' 



z = 2 ULIRSs(Von07) 
SDSS :' 




(b) 



z = ULIRGs(R!eke09) 
SDSS 




'I I I I I' 



(d) 



z = 2.8 ULIRGs(Volionte07) 
SDSS 




10 



10 



5 6 7 8 9 

wavelength(/xm) 

Fig. 2. — The continuum- subtracted comp osite spectrum of the whole SDSS samp le (solid lines) comp ared to those (dotte d lines) of local 
galax ies l lSmith et al.ll2007h , local ULIRGs IIRieke et al.|[2009l). ULIRGs at z =1.7 l IFarrah et al.ll2008r i, ULIRGs at 2=2.8 IValiante et al.l 
l2007t ) and starburst-dominated ULIRGs at z ~ 2 from lYan et alj |(200W) . 



subtracted composite spectrum of the whole SPSS sam- 
ple to those of local galaxies (ISmith et al.1 12007| ). lo- 
cal ULIRGs jlRieke et al.1 f2009h. ULIR Gs at z =1.7 
(iFarrah et aTll2008D , ULIRGs at z=2.8 (jValiante et all 
20071) and starbu rst-dominated ULIRGs at z ~ 2 from 



Yan et all (|2007j ). The composite spectra in the litera- 



ture were used directly while the one of SDSS quasars 
was derived th r ough the arithmetic mea n method (see 
IShi et all 120071 : IVanden Berk et al.l 120011 ) . The contin- 
uum of the SDSS composite spectrum is determined 
by a power law fit to two spectral regions 5.9-6.0 
and 6.7-6.8 /im. For the comparison, the spectra of 
other samples were matched to the SDSS one to have 
the same 6.2 /im aroma tic flux. As s hown in Fig- 
ure EJd), the ULIRGs of IValiante et al.l (|2007h have a 
/(7.7/miPAH)//(6.2/miPAH) ratio consistent with that 
of the SDSS quasars. A slightly higher flux density 
around 7.7 /im for the SDSS quasars relative to lo- 
cal galaxies (Figure ETa))~ local ULIRGs (Figure Efb)) 
and the ULIRGs of lFarrah et~aTl ((2008) (Figure Wc)) 



is most likely due to the presence of the emission line 
[NeV]7.65/im that is absent in star- formation- dominated 
;alaxi es. Only the ULIRGs at z = 2 from lYan et all 
20071) have a /(7.7^mPAH)//(6.2/zmPAH) ratio signif- 
icantly different from the SDSS quasar sample. This 
large deviation is at least partly caused by absorption 
features around 6 /im, such as water ice and hydrocar- 
bons. Therefore, the relative strength of the 6.2 and 7.7 
/xm aromatic features most likely does not change as a 
function of redshift or object type. This implies star- 
formation excitation of the aromatic feature in the SDSS 
quasars. 

In Figure[3j we compare the composite spectrum of the 
whole SDSS sample to those of other quasar samples that 
are selected with different techniques and at different red- 
shifts. For the comparison, all spectra are normalized by 
the mean flux density between 5.8 and 6.1 /im. The 
PG objects are UV/optically-selectcd quasars and thus 
resemble the SDSS quasars. The comparison between 
these two samples should provide clues to the redshift 



9 



E 



1.4 

1.2 

1.0 

0.8 

0.6 
1.4 

1.2 

1.0 

0.8 

0.6 
1.4 

1.2 

1.0 

0.8 

0.6 



(o) 



PG(Shi07) 
SDSS 



(c) 



•+-t- 



(d) 



PG(Shi07) 
SDSS 




X-roy(Sturm06) 
SDSS 



(b) 




PG(Netzer07) 
SDSS 




mm-bright(Lutz08) 
SDSS 



10 



10 



8 9 

wavelength(/xm) 



Fig. 3. — The composite spectrum (solid lines) of the whole SDSS sample compared to those of other quasar samples (dotted lines). 



evolution of the IR properties in UV/optically-selected 
quasars. Figure [Ha) compares the composite spectrum 
of the complete SDSS sample to that of the com plete 
PG quasar sample at z<0.5 from lShi et all (|2007f ). The 
power-law slope of the IR continuum is almost the same. 
The EW of the 6.2 /im aromatic feature is also similar 
in the two samples. In the spectral range of the 7.7 /im 
feature, the SDSS quasars show slightly higher flux den- 
sity, which may be caused by stronger silicate emission, 
IR continuum or 7.7 /im feature emission. As shown in 
Figure OJb) , the su b-sample of the PG quasars used in 
the Q UEST project (|Netzer et al.ll2007t ISchweitzer et all 
2006) shows smaller aromatic feature EWs, which is most 
likely due to their on average higher nuclear luminosity. 
To remove the dependence of the aromatic feature EW 
on the quasar luminosity, we compare the whole SDSS 
sample t o the bright PG sub-sample that is at Mb > 
-23 from lShi et all (|2007| ). The median 5-6 [im contin- 
uum luminosity is 7.5 x 10 10 and 6x 10 10 L for the SDSS 
and the bright PG sub-sample, respectively. As shown 
in Figure [3Jc), the bright PG sub-sample has a 7.7 /im 



feature much weaker than the SDSS quasars while its 6.2 
/im feature is relatively noisy but most likely weaker than 
the SDSS quasars. This indicates that the fraction of the 
aromatic feature flux in the mid-IR evolves with redshift 
at a given nuclear luminosity, which can be caused by 
either higher SFRs in SDSS quasars or the evolution of 
the aromatic feature with redshift. As shown below, this 
evolution is actually dominated by the second effect, i.e., 
stronger aromatic flux for a given SFR at high redshift. 

In Figure [3](d) , we compare the SP SS quasars to the 
mm-bright type 1 quasars at z~2 from lLutz et alj (|2008l ) 
that have mid-IR continuum luminosities two orders of 
magnitude higher than for the SDSS quasars. The IR 
continuum slope and EWs of the two aromatic features 
are similar to those of the SDSS quasars. This again 
suggests that the aromatic feature profile does not change 
with redshift, luminosity or object type. 

The com posite spectrum of ty pe 2 QSOs selected in the 
X-ray from ISturm et all (|2006[ ) is shown in Figure [He). 
This sample spans the redshift range from 0.2 to 1.38 
and has a mean redshift of 0.75. The IR continuum is 
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Fig. 4. — The Lsfir vs - ^6.2umPAH for high-redshift star- 
forming galaxies as listed in Table [2] The solid line and dotted line 
show the linear regression fit to the data and the associated scat- 
ter, respectively. The dot-dashed line s hows the local relat ionship 
based on the star-formation template of Riekc ct al. (2009). 

redder than that of the SDSS type-1 quasars. Although 
no aromatic bands are detected, the spectrum is so noisy 
that even higher EWs than for the SDSS quasars are 
still possible. In addition, the median L^^^n of this X- 
ray-selected sample is about 2.5 xlO 11 L Q , several times 
higher than that of t he SDSS sample. Therefore, in con- 
trast to lSturm et~aU J2006), we argue that the S/N of the 
spectrum is not high enough to rule out the possibility 
that these X-ray-selected type 2 QSOs have significant 
star-formation activity. 

4.2. Star Forming Infrared Luminosity Function of 
Quasar Hosts at z=l 

4.2.1. Conversion factors from the 6.2 [im aromatic flux to 
the total star-forming IR luminosity 

High-rcdshift star-forming galaxies appear to show 
evolution in the relative strengths of the total star- 
forming IR and th e aromatic luminositi es (£sfir/£pah) 
(jRigbv et al.H200a iMurphv et al.ll2009f) . To convert the 
aromatic flux of quasar hosts to the corresponding SFR, 
the ratio (£sfir/£pah) should be defined based on high- 
redshift star-forming galaxies. A heterogeneous sample 
of high-redshift star- forming galaxies at 0.5 < z < 1.5 
was compiled from the literature, as listed in Table O 
The object is named by the corresponding reference fol- 
lowed by the name in the reference. The IRS spectra 
were retrieved from the archive and re-reduced except 
for Teplitz07-2 for which the published spectrum was 
used. Since the aromatic fluxes measured with differ- 
ent m ethods can vary up to several times (jSmith et al.l 
2007), the 6.2 /zm features of these star-forming galax- 
ies were measured through the same method as for the 
quasars. The large equivalent width (> 0.3 /im) of the 6.2 
(im feature in these objects indicates the dominant com- 
ponent of the IR radiation is powered by star formation. 
This is derived from the fact th at the composite H ll-likc 
star-forming galaxy spectra in iSmith et al.l (|2007[ ) have 
EW(6.2/miPAH) of -0.3 /mi. 

To derive the total IR luminosity of these star-forming 



galaxies, their broad-band IR flux densities are fitted 
with the 14 star-forming templates for the IR l u minos - 
ity range of 10 9 75 - 10 13 L Q from iRieke eFall (l2009h . 
As shown in Table H the MIPS 24, 7 and 160 /im 
photometry are used for the lBrand et al.l (|2008h sample, 
while IRS 16 / i m, M IPS 24 and 70 /im are used for the 
IMurphv et all (|2009f) sample. The L S FiR(8-1000/im) is 
measured by integrating the template giving the small- 
est x 2 . For each object, the observed photometry and 
associated uncertainties are Monte-Carlo simulated and 
re- fitted with the star- forming templates. The standard 
deviation of the resulting LsFiR(8-1000/im) is adopted 
as the l-cr error of £sfir due to the observed flux uncer- 
tainties. An additional error in Lsfir due to the scat- 
ter of the template itself is estimated as (Lg C j^ latcI+1 - 
^SFiR atcl ~ 1 )/2, where the template i is the one giving the 
minimum y 2 . The two errors are added quadratica lly to 
give the final error in £sfir- IMurphv et al.l (|2009f ) have 
shown that the ^sfir using photometry including MIPS 
70 or 160 /im flux densities is an unbiased estimate of the 
total IR luminosity based on photometry from the near- 
IR to the sub-mm. For one object (Teplitz07-2) with only 
24 /im data, its 24 /im flux density is converted to the 
total IR luminosity using the redshift-dependent I/24/im- 
-^SFiR re lationship and the l -cr scatter is assumed to be 
0.2 dex (jRieke et al.ll2009l ). IMurphv et all (|2009f ) show 
that, at z ~ 1, the total IR luminosity estimated from 
the MIPS 24 /im flux density does not suffer from any 
systematic shift because the aromatic flux contributes a 
small fraction at rest-frame 12 /im. 

Figure 0] shows the £sfir vs. Lpah for these high- 
rcdshift star-forming galaxies. As a comparison, the lo- 
cal relationship is shown as th e dash-dotted line based on 
the star-forming template of IRieke et alJ (|2009f L There 
is a factor of 2-3 decrease in the isFm/ie^PAH at z ~ 1 
compared to t he local beh avior. We used the IDL code 
linmix_err.pro (Kelly 2007) that employs a Bayesian ap- 
proach to perform linear regression of two variables with 
measured errors. The following relationship is derived 
for quasar hosts at z ~ 1: 

Log(L SF i R ) = (12.16 ±0.09) + 

(1.42 ± 0.29)(Log( J L 6 . 2P A H ) - 9.85) 
±0.25 (1) 

Note that the relationship is obtained between 
Log(LsFm) and L 6 .2PAH - (£6.2PAh)- In this case, the 
slope and intercept of the relationship arc independent of 
each other. The median and 1-a range of the regression 
line are shown as solid and dotted lines, respectively, in 
Figure [U As shown in Table O we do not have star- 
forming galaxies at z ~ 1 with total IR luminosities be- 
low 10 11 Lq. It thus may be invalid to extrapolate this 
relation below 10 11 L®. 

As our objects have LIRG-level SFIR luminosities (Ta- 
ble [1]), the conversion factor isFiR/^PAH is about three 
times lower compared to the local relationship. There- 
fore, our estimate of evolution of star formation in quasar 
hosts is conservative. If we adopt the local relationship 
for our z — 1 quasars, as some works provide tenta- 
tive evidence for no evolution in Lsfir /-k pah ratio (e.g. 
iMagnelli et al.ll2009h . then the evolution from z — 1 to z 
~ would be increased by a factor of about three above 
our estimate. 
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TABLE 2 

High-Redshift Star-Forming Galaxies 



Sources 


z 


-fl6/.im 








Ref. 


^6.2PAH 


EW6.2PAH 








[mjy] 


[mjy] 


[mjy] 


[mjy] 




[10- 14 erg/s/cm 2 ] 


H 


[10"L Q ] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


BrandOS- 70bootes3 


0.986 




1.26±0.1 


35.0±6.3 


145±29 


1 


0.96±0.07 


0.48 


46.39 


Brand08-70bootes4 


0.975 




1.22±0.1 


36.6±5.1 


100±20 


1 


0.84±0.08 


0.40 


41.29 


Brand08-70bootes7 


0.664 




2.18±0.1 


51.8±4.7 


135±27 


1 


1.66±0.08 


0.51 


20.20 


BrandOS- 70bootes9 


0.668 




3.46±0.1 


67.2±3.1 


245±49 


1 


3.29±0.12 


0.79 


29.68 


Murphy09-ID3 


0.63 


0.774±0.006 


1.210±0.005 


11.1±0.53 




2 


1.18±0.03 


0.72 


6.60 


Murphy09-ID8 


0.64 


0.399±0.006 


0.721±0.005 


11.1±0.53 




2 


0.56±0.02 


0.49 


4.93 


Murphy09-ID11 


1.22 


0.993±0.006 


0.446±0.005 


13.2±0.53 




2 


0.76±0.03 


0.73 


53.28 


Murphy09-ID22 


0.64 


0.580±0.006 


0.750±0.005 


5.53±0.53 




2 


0.89±0.03 


0.58 


4.16 


Teplitz07-2 


1.09 




0.133±0.015 






3 


0.16±0.03 


0.70 


3.82 



Note. — References: 1 - IBrand et al.l J2008I) ; 2 - IMurphv et al.l II2009D ; 3 - ITeplitz et al.l J2007D 



4.2.2. Star-forming IR luminosity function of z ~ 1 Type-1 

Quasars 

The SFIR LF of the z ~ 1 type-1 quasar hosts is 
constructed by counting the number of objects in dif- 
ferent luminosity bins and redshift bins (ALAz). For 
this method, in the lowest luminosity bin that intersects 
the flux limit at a given redshift, only the part of the 
ALAz space above the flux limit can be filled by ob- 
jects. To account for this partly filled bin, we ha ve used 
a revised 1/Vm ax method (|Page &; Carreral [20001 ) to de- 
rive the luminosity function. The inco mpleteness func- 
tion fo r the quasar selection is given by Richard s et alj 
(12006( 1 . The luminosity bins are defined to have roughly 
the same number of objects and range from the observed 
lowest luminosity to the highest one. 

Figure [5] shows the SFIR LF of z ~ 1 type-1 quasar 
hosts compared to th e LF o f contemporary field galaxies 
from iLe Floc'h et all (|2005l ) and to that of the z=0.25 
PG quasars. The SFIR LF of the PG quas ars has been 
updat ed by using the £sfir/£pah ratio of iRieke et al.1 
(2009), but there is only a small change. The error 
bar shows the Poisson uncertainty. The uncertainty in 
the SFIR luminosity can smooth the LF by scattering 
more objects out of the luminosity bin originally with a 
larger number of objects. To quantify this effect, a set 
of 1000 SFIR LFs were Monte-Carlo simulated using the 
observed Lsfir and associated uncertainties. The rela- 
tive change in the number density for each luminosity 
bin is then quantified. The grey area then shows the 1-cr 
range of the SFIR LF after correcting for the effects due 
t o the unce r tainty of the SFIR luminosity. 

IShi et all (|2007l ) carried out a series of simulations to 
demonstrate that the PG quasars have a flat SFIR LF 
compared to that of field galaxies. As shown in the figure, 
such a trend most likely exists at z ~ 1. Within our 
limited range of luminosity, the difference between the z 
~ 0.25 and z ~ 1 SFIR LFs can be described as pure 
density evolution. 

4.2.3. Cosmic Evolution of the Comoving SFIR Luminosity 

Density And Average SFR in Quasar Hosts 

With the SFIR LF at z ~ 0.25 and z ~ 1, we now 
can quantify the cosmic evolution of the comoving SFIR 
luminosity density in type-1 quasar hosts. Due to the 
lack of constraints on the slope at low-luminosity, two 
types of estimation have been explored. First, we can es- 
timate the comoving SFIR density by simply integrating 



the available data points, which gives an evolution factor 
of 10±4. Second, we can fit the SFIR LF with a series 
of Schechter functions whose slopes (a) are Monte-Carlo 
simulated. For the PG sample, the slopes are simulated 
using the number densities and associate uncertainties in 
two low luminosity bins. For the SDSS sample, as the 
slopes given in a similar way have large uncertainties, we 
attempt to use the objects with undetected aromatic fea- 
tures to constrain the slope a. Basically, for a range of 
assumed slopes, we fitted data points to obtain a series 
of analytic Schechter LFs. For each LF, a SFIR lumi- 
nosity is then randomly assigned to a SDSS quasar with 
the relative probability following that of the LF. The de- 
tection rate of the aromatic feature is then calculated 
by comparing the simulated SFIR luminosity to the ob- 
served 3-cr uncertainty of the SFIR luminosity. Figure [6] 
shows the simulated aromatic detection rates (points) 
compared to the observed one (lines) . The derived slope 
is a = — 0.59±0.25. By Monte-Carlo simulating this 
slope, the observed LF is fitted with a series of Schechter 
functions, which is shown as a solid line in Figure [5] and 
listed in Tabled We obtained Logi* = 11.64±0.19. The 
comoving SFIR luminosity density is then obtained by in- 
tegrating the fitted Schechter functions and the result is 
shown in the Figure [7£a). In this method, the derived in- 
crease of the comoving SFIR luminosity density in quasar 
hosts from z—0.25 to z=0.9 is 15^g (1-c)- In either case, 
the SFIR energy density in type 1 quasar hosts shows a 
dramatic evolution, much larger than the behavior of the 
general field galaxies (which increase by a factor of 5^J 7 ; 
ILe Floc'h et ani27)05l[Perez-Gonzalez et al.ll2005h . 

The integration of the LF is carried out in the luminos- 
ity range [10 10 , 10 12 ] L©. If we vary the upper- or lower- 
limits by a factor of 10, the derived evolution changes 
little. If we had used the local conversion of the aro- 
matic luminosity to the SFR, the evolution of the quasar 
host SFRs would roughly triple. 

We now consider how selection biases might influ- 
ence this result. We quantified the evolution of the 
comoving SFIR luminosity density in quasar hosts by 
measuring the SFR in PG quasars and SDSS quasars, 
two samples selected through different UV/optical cri- 
teria. The SDSS quasars were selected on the basi s 
of multiple color-color diagrams (Ric hards et alj l2001h ; 
nonetheless, as with the PG sample, UV excess is the 
dominant characteristic. The incompleteness of the PG 
quasar selection (U-B < -0.46 and B < 16.16) has long 
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Fig. 5. — The star-forming IR luminosity function (SFIR LF) of z ~ 1 quasar hosts (filled circles) compared to that of low-redshift PG 
quasars (open diamonds). The grey area shows the 1-cr range of the SFIR LF after correcting the effect due to the uncertainty of the SFIR 
luminosity (see text for details). The SFIR LFs of field galaxies at both low and high redshifts are also shown for comparison. 



TABLE 3 

Best-fitting parameters to star-forming IR LF of SDSS and PG 



QUASARS 


Object Log( 


^>*[Mpc- 3 LogL- 1 ]) 


Log(L*[L Q ]) 


a 


SDSS 


-6.30±0.16 


11.64± 0.19 


-0.59±0.25 


PG(M S <-23) 


-7.86±0.39 


12.04±0.34 


-1.10±0.19 



Note. — The slope for the SDSS sample is derived by producing the 
observed PAH detection rate (sec Figure. [6j, while the one for the PG sample 
is calculated using the number densities and associate uncertainties in two 
low luminosity bins. By Monte-Carlo simulating the derived slope, the cf>* 
and L* is obtained by fitting Schcchter functions to the observed LFs. All 
uncertainties are given at 1-a level. 
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Fig. 6. — The triangle shows the predicted detection rate of the 
aromatic feature of the SDSS sample for a range of the assumed 
Schechter function slope. The solid line is the observed detection 
rate while the associated uncertainty is shown as the dotted line. 

been known dGoldschmidt et all Il992t IWisotzki et all 
I2000t iMickaelian et all 120011 1 Jester et all l2005fl . Com- 
parisons of PG quasars to other quasar surveys indi- 
cate a completen ess of ~50-100% (I Wisotzki et "all 120001: 
IMickaelian et alj I200H I Jester et al J 120051 ). Theincom- 
pleteness is most likely independent of the quasar op- 
tical and radio pro perties within the PG quasar selec- 
tion criterion itself l|Jester e t al. 2005), i.e., PG quasars 
are representative of quasars that are bright at B-band. 
Therefore, the evolution of the comoving SFIR luminos- 
ity density in type 1 quasar hosts is over-estimated at 
most by a factor of 2. In addition, our z ~ 1 quasars 
(Mi < -24.7; M B <~ - 24.0) are brighter than our low- 
rcdshift PG quasars by one magnitude (Mb < -23). Our 
study of PG quasars showed that brighter PG quasars 
have a flatter SFIR LF, implying on average a stronger 
SFR in a brighter PG quasar host. However, the inte- 
gral SFIR density is actually lower by a factor of 2 for PG 
quasars at Mb < -23 compared to the value for Mb < 
-22, simply due to lower number densities. The evolu- 
tion of the comoving SFIR density in quasar hosts may 
be even larger if we compare quasars at the same bright- 
ness cut. Therefore, the two factors (in-completeness of 
the PG sample and un-matched brightness limit of two 
quasar samples) tend to cancel each other out. We con- 
clude that the evolution of the comoving SFIR luminosity 
density in type-1 quasar hosts is indeed much larger than 
that in field galaxies. 

We can use the SFIR LF of quasar hosts to derive the 
average SFR per object: 

/mD v /o°° SFR Lsfir $sfircLL S fir , . 

(blK qS o) - p5o "77- > i 2 ) 

J j limit ^qso^-^qso 

where $sfir is the SFIR LF of quasar hosts 
and $ qso is the LF of the quasar luminosity at 
the wavelength where the quasar sample is se- 




0.4 0.6 
redshift 



Fig. 7. — The cosmic evolution of the comoving SFIR luminosity 
density (a), the average SFR (b) and the average SFR/M acc (c) 
in type 1 quasar hosts . The SFIR luminosity is converted to the 
SFR u sing Ricke et al. (2009) relations that deviate from those of 
Kcnnicutt (1998) by a factor of ~2 as a result of a different initial 
mass function. 



lected. The result is shown in Figure \T[b), where 
isFiR is converted to the SFR using the rela- 
tion SFR(M /yr)=6.37xlO- 12 (LTiR/i0) l n at L Tm > 
10 11 L and SFR (M^/yr)=2 31xl Q- 11 aTTR/^) 106 at 
Ltir < 10 11 L (IRieke et al.H2009D. T hese relationships 
deviate from those of iKennicuttl ( 19981 ) mainly due to a 
different initial mass function. As shown in Figure [TJ^b), 
the average SFR in the quasar host is LIRG-level (and 
remains in the LIRG range with the local isFiR/^PAH 
applied at z ~1). The evolution of this average SFR is 
very insignificant, and remains small even with applica- 
tion of the local £sfir/£pah a.t £ ~ 1. As we noted 
above, our SDSS quasar sample is limited at a mag- 
nitude (Mb <-24) brighter than the PG quasar limit 
(Mb <-23). Brighter quasars have on average higher 
SFRs (|Shi et al.ll2007l ). PG quasars at M B < -23 have 
SFRs on average a factor of 1.3 higher than PG quasars 
at Mb < -22. Therefore, the average SFR in the optical 
type-1 quasar is most likely nearly constant with red- 
shift. This result supports our conclusion that most of 
the host galaxy SFIR LF evolution is in density, not in lu- 
minosity. Here we simply assume that the same SFR- Mr 
trend holds at z=0 and z=l and that the trend between 
Mr=-22 to -23 holds up to Mr=-24, which is most likely 
true based on the correlation between SFR and quasar 
nuclear luminosity as discussed in § 14.31 Note that the 
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derived quantities in Figure [7] do not depend on the in- 
tegration limit. For example, the change is only about 
0.1 dex using the observed range of the SDSS sample 
([10 n ,10 12 ]). 

4.3. The SFIR Luminosity As A Function of Quasar 
Nuclear Luminosity 

The SFIR LF of the type-1 quasar hosts gives a global 
sense of their star-forming activity. The enhanced SFR 
in the quasar hosts as implied by the SFIR LF indicates 
that further investigation of the star formation activity 
in individual objects is critical to understand the inter- 
play bet ween star forma tion and nuclear activity. As 
shown in lShi et alj (|2007| ). the brighter PG quasars have 
on average higher SFRs as measured from their SFIR LF. 
We now search for direct correlations between the SFIR 
luminosity and quasar nuclear luminosity. 

Figure [8] shows the result for the whole PG sample of 
90 objects (triangles) and the whole SDSS sample of 57 
objects (stars). As there is no significant redshift evolu- 
tion of the SFRs in the quasar hosts (§ I4.2.3[) . the two 
quasar samples are plotted together. The 5-6 /jm IR 
continuum luminosity {L^gl^A is used as a tracer of nu- 
clear luminosity. For clarity, the uncertainty of the SFIR 
luminosity (Xsfir) is not shown in the figure. A signif- 
icant number (50%) of objects only have upper-limits. 
To demonstrate there is a correlation between the SFIR 
luminosity (Lsfir) and nuclear luminosity, we stacked 
spectra of several sub-samples of quasars without aro- 
matic detections defined within certain Lf^g luminos- 
ity bins. The result is shown in Figure [31 Each compos- 
ite spectrum has detected aromatic features. The 6.2 /im. 
feature is measured as in this paper while the 7.7 and 11.3 
/zm features are measured following lShi et all (2007). For 
each sub-sample, the mean L^R m is used and the cor- 
responding SFIR luminosity is measured based on the 
composite spectrum. The result is shown as large sym- 
bols in Figure [5] The m range of each sub-sample 
is used as the uncertainty of the and the corre- 
sponding -Lsfir range is used as the uncertainty of the 

^SFIR- 

To quantify the correlation, linear regression fits were 
carried out using linmix_err.pro that accounts for the 
measured errors and upper- limits. We performed fits to 
three data sets: all individual PG+SDSS data points, all 
detected + stacked data points and all PG objects. The 
best-fit lines are shown in Figure [8] and listed in Table [H 
As shown in the figure, the three fits give roughly similar 
results and have correlation coefficients of 0.6-0.8. Such 
consistent results confirm the existence of the correlation 
and the lack of systematic difference between z = 1 SDSS 
and low-redshift PG samples. 

To further demonstrate that the correlation is not 
caused by the aromatic feature non-detections in about 
half the objects, we carried out a Monte-Carlo simulation 
for the PG sample, which spans a large range of SFIR 
luminosity. Basically, a SFIR luminosity is assigned ran- 
domly to a PG quasar with the relative probability fol- 
lowing the SFIR LF of the whole PG sample. In such a 
simulation, the SFIR luminosity is assumed not to cor- 
relate with the nuclear luminosity. For each set of sim- 
ulations, a sub-sample of objects with aromatic detec- 
tions can be defined through comparing the simulated 



^SFiR and the observed uncertainties or upper- limits. 
The Spearman's correlation coefficient can be measured 
for the Log(LsFiR)— Log(L^ u ^ tm ) of such a sub-sample. 
After 10000 simulations, the probability of the simulated 
Spearman's correlation coefficients being higher than the 
observed one is only 1%. 

Both of these tests imply that the SFIR luminosity 
most likely correlates with the nuclear luminosity in type- 
1 quasar hosts. The comparison of this correlation to 
those in theoretical models is presented in § 15.31 



5. DISCUSSION 



5.1. 



Comparison with Other Studies of AGNs at 
High- Redshift 



Intense star formation in high-redshift quasars has 
been observed in some individua l quasars as indica ted by 
rest-frame far- IR emission (e.g. IWang et all [20081. aro- 
matic features dLutz et al.ll2008fh molecular gas em ission 
(for a review, see ISolomon fc Vanden Boutf [2005) and 
UV emission ( Akivama |2005|) . However, these studies 
suffer from low number statistics, selection bias toward 
high SFRs and large uncert ainties in the derived SFR. 
iHatziminaoglou et al.l (|2008l ) have studied IR SEDs of 
SDSS quasars in the SWIRE field through fitting torus 
models and star- formation templates. They derived on 
average star-formation contributions of ^35% to the to- 
tal IR luminosity for quasars at z =0.8-1.0, correspond- 
ing to a SFIR luminosity of ~ 10 12 L Q . However, this 
high SFR is only for objects detected at MIPS 70 and/or 
160 /im, indicating a bias toward high SFRs in their 
sample. Our study of a complete sample of 57 SDSS 
quasars demonstrates conclusively that type-1 quasars 
reside in host galaxies with intense star formation. As 
the nucleus is brighter due to selection effects and the 
host galaxy is apparently smaller and fainter due to 
surface-brightness dimming, the host morphologies are 
difficult to c onstrain unambiguously (for a review, see 
lDaviesll2008D but appear to be consistent with elliptical 
morphologies (iKukula et al.ll200il; iHvvonen et al.l 120071: 
Falomo et alJl2004t iKotilainen et al.ll2007t iFalomo et all 
2008t I Amnions et al.ll2009[ ). We will discuss this appar- 



ent inconsistency between early-type host morphology 
and intense star formation in § 15.2.11 

Similar to their low-redshift counterparts, IR-selected 
type 2 quasars l argely show distur bed morphologies 
(|Lacv et al.l 120071 : lUrrutia et al.l 120081 ) and have intense 
star formation act i vity as probed by the aromatic feature 
(|Lacv et al.l l2007t iHernan-Caballero et alj I2009D . High 
SFRs are also found for optically- selected type 2 quasars , 
based on the aromatic feature (Za kamska et al.l 12008). 
High-z ULIRGs with embedded AGNs (possibly related 
to type 2 quasars) have also been f ound to be associ- 
ated with intense star formation (e.g. iHouck et al.ll2005 



Sajina et al.l l2007 



I Weedman et al.l[200H Yan et alJl2007t 
IWatabe et al.ll2009f ). 

A correlation between the mean stellar population 
age and the AGN lu minosity is found for h igh-redshift 
X-ray-selected AGNs ([Amnions et al.ll2009f ). Although 
some studies do not confirm strong on-going SFRs in 
X-ray-sel ected AGNs with a range of the X-ray lu- 
minosity dSturm et al.|l2006t [Alonso-Herrero et al.l 120081 : 
ISchawinski et al.1 |200 ff), a more thorough investigation 
of star formation in a relatively large sample with more 



15 



TABLE 4 

Best-fitting parameters to the correlation i/SFiR-£ NUC 



data sets 




Best-fitting Formula 








correlation coefficient 


All PG + SDSS 
All detected plus stacked points 
All PG 


Log(LsFm) 

Log(LsFIR) 

Log(LsFiR) 


= (10.98 ± 0.05) + (0.97 ± 0.08)(Lo f 
= (11.00 ± 0.07) + (1.14 ± 0.11)(Lof 
= (10.54 ± 0.12) + (0.86 ± 0.15)(Lof 


t(I>6 


-6/im) 
-6/im) 
-6fim ) 


- 10.58) ±0.52 

- 10.42) ±0.56 

- 10.37) ±0.75 


0.76±0.05 
0.81±0.05 
0.64±0.11 



Note. — The correlations arc derived between Log(LgFiR) and LogL 
and slope are independent. 
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Fig. 8. — The correlation between the star-forming IR luminosity and the nuclear luminosity at 5-6 /an. The open triangles and stars 
are for individual PG and SDSS quasars, respectively. The large filled symbols show the result based on the composite spectra of PAH- 
undetected objects within 5-6 fim luminosity ranges as indicated by their error-bars on the X-axis (see Figure [5J. Three linear regression 
fits are listed in Table l4l 



plausible probes should be made before concluding that 
the X-ray-selected AGNs are in a different evolutionary 
stage. For example, a recent study of X-ray-selected 
AGNs in the COSMOS field indicates enhanced SFRs in 
their host gala xies relative to field ga laxies with similar 
stellar masses (jSilverman et al.l f2009). It also suggested 
that a selection bias may account for the X-ray-selected 
AGNs residing in the green valley region of the color- 
magnitude plot. 

To summarize, AGNs at high redshift (z > 0.5), in- 
cluding high-luminosity type-l/type-2 quasars and rel- 
atively low-luminosity X-ray-selected AGNs, are experi- 
encing intense star-forming activity in their host galaxies. 



5.2. Quasar Host Galaxies As A SFR-Enhanced Subset 
of Field Galaxies 

We have measured the aromatic-based SFR for a com- 
plete sample of SDSS type 1 quasars at z ~ 1. The 
derived SFIR LF is flatter than that of z ~ 1 field 
galaxies, implying enhanced SFRs in the quasar hosts. 
Combining this with studies of IR - and optically-selected 
type 2 quasars (lLacv et alJl2007t IZakamska et alJl2008t 
iHernan-Caballero et al.l l2009f ). we conclude that the 
quasar host galaxy population is a SFR-enhanced subset 
of field galaxies. This has been known to be true for low- 
rcdshift quasars as shown in our study of local optically- 
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Fig. 9. — The composite spectra of PG and SDSS quasars without 
PAH-detection within given 5-6 (im luminosity ranges. 

selected PG, IR-selected 2MASS and radio-selected 3CR 
quasars (|Shi et al.ll2007ft . 

As discussed in the introduction, there are many differ- 
ent ways to probe the stellar population age. To interpret 
results with different probes, one must keep in mind that 
each probe is only sensitive to stellar populations with 
certain ranges of ages and suffers from different limita- 
tions. Therefore, apparently different results about the 
stellar population age of quasar hosts by different studies 
are not necessarily inconsistent with each other, but in- 
stead should be integrated to achieve the complete view 
of the stellar population that is critical to understand the 
BH/galaxy co-evolution. In the following, we reconcile 
our result of enhanced SFRs with different results in the 
literature. 

5.2.1. Comparisons with Host Galaxy Morphology Studies 

The IR- and optically-selected type 2 quasar 
hosts show high merger (>50%) fractions at both 
low (ICanalizo fc Stockton] 120011: IHutchings et all 120031: 



Marble et al. 2003) and high redshift (lLacy et al.1 



20071 lUrrutia et al.l 



2008; iHernan- Caballc ro et al.l l2009) , 



higher th an field galaxies of similar infrared luminosity 
(~30%; iShi et al.l 120091 iSobral et al.l [2009h . The en- 
hanced SFR in type 2 quasars is thus consistent with 
their host morphologies. 

For optically/UV-selected type 1 quasar hosts, 
regular early-type morphologies dominate, with a 
fraction of >5 0% for quasars at M B < -23 at low red- 
s hift (z<0.5) (IHutchings et al.l 119841: iSmith et al.l 11986' 



shift, the host galaxi es appear to be cons i stent with ellip- 
tical m orphologies (iKukula et all 1200 It iHvvonen et al.l 
20071: iFalomo et al l 12004 iKotilainen et alj 120071: 
Falomo et al.l l2008t lAmmons et al.l I2009D . despite the 
difficulties in constraining the types unambiguously. 
The merger fraction is <30%, not significantly different 
from field galaxies. Therefore, the enhanced SFR in the 
type 1 quasar hosts is not a natural result of the fact 
that the host galaxy is massive, given the dominance of 
regular early-type morphologies. 

To reconcile the apparent inconsistency between 
elliptical-dominated morphologies and enhanced SFRs, 
we propose that star formation in type 1 quasar hosts 
occurs in the circum-nuclear (^0.1-lkpc) region. In 
this case, any star-forming signature would be missed 
by spatially-resolved image studies that always ex- 
clude the central region as a result of subtracting the 
diffracted nuclear light. The large extinction in the cir- 
cumnuclear region can also hide star formation traced 

by extinction-sens itive probes (e .g. UV emission & 

[OII1 A3727 lines) (|Akivamal 12003 : iHol l200l iKim et all 
2006). Current examples of direct detections of in- 
tense circum-nuclear starbursts in type-1 quasars ar e 
limited to a few nearby objects (|Cresci et al.l [2004). 
High-resolution molecular gas mapping has revealed 
a large gas concentrat ion in the circum-nuclear re- 
gion for tens of quasars (ISolomon &, Vanden Boutll2005l: 
Maiolino et~aT1 120071 : iRiechers et all 120091 : IWalter et all 
2009( 1. As the low-luminosity counterparts of quasars, 



nuclear starbursts (Heckman et aI][T997: Gu et al.l 2001: 


[Gonzalez Delgado et al.ll2001l: 


Cid Fcrnandcs et al. 


2004: 


bavies et al.ll2007t iRiffel et al. 


2009), although enhanced 



star format ion also occurs in th e spiral disks of Seyfert 
2 galaxies (|Maiolino et al.l lT995). Quasars, with higher 
luminosity and thus larger mass inflow, probably harbor 
more intense circum-nuclear starbursts. As shown in Fig- 
ure E£b), quasar hosts have on average LIRG- level SFRs. 
In normal local galaxies, such intense star formation is 
normally ach ieved only in th e circumnuclear region (for 
a review, see iKennicuttl Il998h . 

5.2.2. Comparisons with the Study of Stellar Populations 

Studies of stellar populations through broad-band 
SEDs and optical spectra have discovered intermediate- 
age stars or bluer colors in AGN hosts compared 
to their i nactive galaxy co unter par ts at both low 
redshift dRonnback et al.l 119961 : iBrotherton et "ail 

Uahnke fc Wisotzkil 
" iCanalizo et al.l | 2006| : 
Jahnke et all [2007 



1999: 
2003; 



Kauff mann et alj 120031: 



.Jahnke et al.l 
Vanden Berk et ahl 



2001 
I2006T 



McLcod & Rickc 1995; Bahcall ct al. 19971: iDunlop et all 
20031: iFlovd et al.ll2004l : lGuvon et alJhoOffl . At high red- 



Schawinski et al.ll2009f ) and high redshift ([Amnions et al.l 
20091) . In the following, we show that the spectral fea- 
tures employed by current studies of stellar populations 
are not sensitive to on-going (<0.1 Gyr) starburst 
activity. 

The continuum shape is useful to detect the pres- 
ence of the young stellar population only if the spec- 
trum extends shorter than 4000A or the broad-band 
SED includes U-band. At low redshift, it is diffi- 
cult to achieve good coverage of the spectrum/ SED 
at such short wavelengths. Thus our result of en- 
hanced SFRs does not contradict studies based on 
the continuum shape at A > 4000 A that detect 



interm ediate-age stellar populat io ns (|Ronnback et alj 
19961: IJahnke fc , Wisotzkl 12001 IJahnke et all l20f)l 
Schawinski et alj 120091 : I Amnions et alj |2009| )" On the 
other hand, at A < 4000 A the contrast between nuclear- 
and host-light in type 1 AGNs becomes much higher 
compared to the value at redder wavelengths. The sub- 
traction of the nuclear light always introduces large er- 
rors on the resulting host brightness. Even in type 2 
AGNs, the scattered nucle ar UV light is sometim es im- 
portant or even dominates (Zakamska ct al. 2006). More 
importantly, obscuration can hide blue host-galaxy light 
if star formation occurs in the circum-nuclear region as 
discussed above. Therefore, results based on spectra 
extending shorter than < 4000 A may still underesti- 
mate sig nificantly the level o f on-going star- formation 
activity (IJahnke et all l2007t IVanden Berk et alJ 120061 : 
iChen et all 120091). 

iKauffmann et al.l (|2003l ) employed the 4000 A break 
and H<5 absorption to study stellar populations and found 
intermediate-age stars in hosts with early-type morpholo- 
gies. The 4000 A break is caused by metal line absorp- 
tions. In hot stars with age < 0.1 Gyr (see Fig. 2 of 
IKauffmann et all (|2003bh ). the 4000 A break is small 
(< 1 .15) and insensitive to the stellar age. Moreover, the 
4000 A break is a luminosity-weighted mean-age indi- 
cator and thus does not have the temporal resolution 
to separate young stars from dominant old stars in a 
massive host galaxy. H5 absorption arises in late-B to 
early-F stars and i s thus only prominent .1-1 Gyr af- 
ter star formation (|Kauffmann et al"]l2003bl ). Therefore, 
the two probes are not sensitive to ongoing star forma- 
tion activity. In addition, stellar-age-dependent obscu- 
ration may also completely hide the effect of massive 
OB stars on these two features. In fact, such selec- 
tive obscuration may be the reason for the presence of 
some LIRGs whose optical spec tra show strong Ht5 ab - 
sorption but weak 0[II] emission ( Poggia nti fc Wull2000f ). 
Studies based on the whole SPSS spectrum suffer from 
similar problems (IJahnke et al.ll2007tlVanden Berk et al.l 
l27)M I Wild eXall 120071 : IChen et al.ll2009h . as the contin- 
uum shape may be affected by extinction and nuclear 
light contamination while all the available stellar features 
(Balmer break, high order Balmer absorption line, 4000A 
break and Ca II (H&K)) are not sensitive to the current 
star formation whose characteristic features (mainly neb- 
ular emission lines) are contaminated severely by the nu- 
clear radiation in both type 1 and type 2 objects. 

In a summary, by integrating different studies, we now 
have a more complete census of the stellar population in 
type-1 quasar hosts, which appear to be massive galax- 
ies presumably dominated by old stars, but that harbor 
a significant fraction (10%) of intermediate-age stellar 
populations, and are experiencing intense circum-nuclear 
star formation. 

5.3. Implications For the BH-Bulge Correlation 

The direct measurements of SFRs in quasar hosts can 
provide new insights into the mechanism that shapes 
the A/bh-c relationship. Figure [7(c) shows the redshift 
evolution for the ratio of the SFR in quasar hosts to 
the black hole growth rate, obtained by dividing the 
integral SFIR luminosity function by the integral BH 
growth rate function. The SFIR luminosity is con- 
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Log(dM BH /dt[M /yr]) 

Fig. 10. — The ratio of the SFR and the BH accretion rate vs. 
the BH accretion rate. The open circles show the results for our 
2 ~1 quasars and the filled circle is the mean value derived from the 
luminosity function. Curves show these two quantities along the 
m erging process i n two m odels, where the solid line is for the model 
in lGranato et aLl (120041) a nd the dot-dashed line is for the model 
in lDi Matteo et a l. (2005). The thick part of the curve indicates 
the stage after the peak BH accretion rate (presumably the type-1 
quasar phase). 
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Fig. 11. — The correlation between the star-forming IR, luminos- 
ity and the nuclear luminosity at 5-6 (im for the stacked spectra of 
all PG and SDSS quasars within certain 5-6 fim luminosity ranges 
as indicated by their error-bars o n the X-axis. The s e are compared 
to the predictions of the model of Di Matteo ct al. (2005) and that 
of lGranato et aLl (120041 ) shown as dashed and dotted lines, respec- 
tively. 
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verted to the SFR using the iRieke et all (|2009f ) rela- 
tions as shown above. The BH growth rate function 
is obtained from the quasar £?-band luminosity function 
using L b oi = e/(l - e)Mbhc 2 = c B L B , where L hoi is 
the bolometric luminosity, e is the mass to energy con- 
version efficiency, Mbh is the BH growth rate and cb 
is the bolometric correction for the B-band luminosity 
Lb- Figure[7Jc) shows the result for e=0.1 and cb=H-8, 
wh ere the quasar i?-band lu mino sity functions are give n 
bv ISchmidt fc Greenl (|1983h and iRichards et all (j2006f h 
As shown in the figure, the SFR/A/bh in quasar hosts 
is almost constant with redshift and has a value of ~10, 
much lower than that expected from the Mbh-c relation- 
ship. This result rules out the simple model that the BH 
accretion and star formation evolve with a fixed ratio as 
given by the Mbh-c relation at any time for individual 
objects. This implies that Mbh/M* for type 1 quasars 
deviates from that implied by the Mbh-c relation. The 
factor of deviation of Mbh/M* can be given roughly as 
(l+AA/ B H/^BH)=exp(?7 4 4 ^ 1 B ° 7yr ), where r] is the frac- 
tional Eddington accretion rate ry=M acc /Af ™ D and t qso 
is the duration of the type-1 quasar phase in years. For 
r\ = (0.1-1) and £ qso =10 8 yr, the deviation is a factor of 
1-10. 

Numerical simulations have invoked galaxy mergers 
and quasar feed back to successfully e xplain the Mrh- 
a re lationship (|Granato et alj l2004t iDi Matteo et al.l 
I2005T ). Different models predict different time evolution 
of the SFR and BH accretion rate along the merging 
process. Our direct measurements of SFRs in type 1 
quasars can provide constrains on these models. Our 
z ~1 quasars hav e a median BH mass of 10 8 - 9±0 - 3 M Q 
ijShen et al.l(200l. In Figure [TOl we compare our result 
to the model of ldranato et a l. ( 2004) for their most mas- 



sive d ark matter halo and the model of IDi Matteo et all 
(2005) for their most massive galaxy. For each model, 
the thicker part of the curve indicates the stage after the 
peak BH accretion rate, presumably the type-1 quasar 
phase. The open circles show the results for our z ~1 
quasars and the filled circle is the mean value derived 
from the luminosity function. As shown in Figure 1101 
individual quasa rs span a range of Mb h and SFR/A/bh- 
For the model of IGr anato et all (|2004l ). the full range of 
the SFR/A/bh resides within their type 1 quasar phase. 
However, the SFR in the model is more than an order of 
magn itude higher than our average value (the filled cir- 
cle). IDi Matteo et al.l (|2005f) produced a model with the 
mean Mbh and SFR/A/bh comp arable to our observed 
result. The much higher SFR in IGranato et~aTl (|2004D 
may be related to their high star formation efficiency 
which is determined by th e free-fall or coo li ng tim e in 
their models, while in the IDi Matteo et alj (|2005f) star 
formation occurs in a self-regulated mode. 

In Figure [TT1 we further compare our Lsfir - ^{?-6u m 
relationship to predictions of these two models. As mod- 
els only produce the average case, we here use the £sfir 
- Lf^Q^ m relation determined from the composite spec- 
tra of all PG and SDSS quasars within certain luminosity 
ranges. The SFRs of the tw o models are conver ted to the 
SFIR luminosity using the IRieke et all (|2009f) relation- 
ship. The BH accretion rate is first converted to the bolo- 
metric luminosity by assuming Lboi = — e)Mbhc 2 



and e=0.1. The resulting bolometric luminosity is con- 
verted to the based on the Lbol-^s'-eum relation- 
ship of our SPSS samp le whose bolometric luminosity is 
from lShen et al.l (|2008f ). Only the type-1 stage (i.e., after 
the peak BH accretion r ate) of the two models is shown in 
FigureEJ The model of lDi Matteo et al.l <|2005l ) is gener- 
ally consistent wit h our observed relatio nship with a bit 
lower slope, while IGranato et alj (|2004h predict a much 
steeper slope. This implies that t he quasar feedback ma y 
be overestimated signifi cantly in | Granato et a .1 (2004). 
The feedback energy in lDi Matteo et aIT "( 2005 ) is fixed 
to be 0.5 % of the accreted rest mass energy while this 
fraction in IGranato et a l. (2004) is several times higher. 



6. CONCLUSIONS 

To constrain the cosmic evolution of star formation in 
quasar host galaxies, we present Spitzer IRS observations 
of a complete SDSS sample of 57 type-1 quasars at z ~1. 
The main conclusions are: 

(1) About half of the sample has aromatic features de- 
tected at 6.2 and/or 7.7 fim. The composite profile of 
these two features is similar to those of normal galaxies, 
ULIRGs and AGNs at both low and high redshift, im- 
plying star-formation excitation of the aromatic features 
at z ~1. 

(2) Based on the aromatic-to-SFR ratio of star-forming 
galaxies at z ~1, we have constructed the star- forming IR 
(SFIR) luminosity function (LF) of z ~1 type-1 quasars. 
Similar to low- redshift PG quasars, these z <~1 quasars 
show a flatter SFIR LF than z ~1 field galaxies, imply- 
ing the quasar host galaxy population has on average 
higher SFRs. Based on the measured SFIR LF, individ- 
ual quasar hosts are shown to be experiencing on average 
LIRG-lcvcl SFRs, which most likely occur in the circum- 
nuclear region. 

(3) By comparing with similar measurements of low- 
redshift PG quasars, we found that the comoving SFIR 
energy density in type-1 quasar hosts shows much faster 
evolution than that in field galaxies, while the average 
SFR and the average SFR/M acc ratio are almost con- 
stant with redshift. 

(4) For individual objects, we have found a correla- 
tion between the aromatic-based SFR and the nuclear 
luminosity. By comparing this result to predictions 
of different models, we have found that the model in 
IGranato et al.l (|2004h may over-estimat e the quasar feed- 
back significantly while the model in IDi Matteo et al.l 
(2005) produces a roughly consistent result. 

(5) By combining different studies in the literature, we 
now have a more complete view of the stellar population 
in type-1 quasar hosts, which reside in massive galax- 
ies with dominant old stars, harbor a significant fraction 
(~10%) of intermediate-age stellar populations and are 
experiencing intense star formation. 
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